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Abstract 

We study the scaling properties of the clusters grown by the Wolff algorithm 
on seven different Sierpinski-type fractals of Hausdorff dimension 1 < < 3 in 
the framework of the Ising model. The mean absolute value of the surface energy 
of Wolff cluster follows a power law with respect to the lattice size. Moreover, 
we investigate the probability density distribution of the surface energy of Wolff 
cluster and are able to establish a new scaling relation. It enables us to introduce 
a new exponent associated to the surface energy of Wolff cluster. Finally, this new 
exponent is linked to a dynamical exponent via an inequality. 
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In 1966, Kadanoff proposed an intuitive renormalization picture IJ to explain the 
Widom's phenomenological hypothesis This hypothesis states a homogeneous trans- 
formation of the singular part of the free energy per spin under a change of length unit 
from 1 to 6 in the vicinity of the critical temperature and can be written as 

f{t,h) = h-^f{hy%hy-h), (1) 

where d is the integer space dimension of a translationally invariant system; the reduced 
temperature t = T/Tc — 1 and the external field h are supposed to be small; yt and yt are 
two exponents associated with the two relevant directions of the renormalization flows. 
One notices that a translationally invariant system of dimension d is auto-similar and can 
be considered as a particular case of fractal. In the case of a general fractal system, the 
translational symmetry is broken; however, a generalization of the space dimension can 
be done by introducing the Hausdorff dimension df. The number of spins on such system 
can be written as = L'^f where L is the network size. Thus, the Widom's hypothesis 
can be generalized by replacing the factor b~'^ in Eq.(Q) by to describe the decrease 
of the effective number of spins during the change of the length scale. Since fractals are 
constructed by iteration of a generating cell, b cannot take any value; it must be chosen so 
that the fractal structure remains invariant under the change of length unit. The critical 
behaviors of spin models on fractals whose Hausdorff dimension lies between 1 and 3 have 
been studied intensively by Monceau and Hsiao and co-workers [21 IH El El and Carmona 
et al. [71, by performing Monte Carlo simulations. The results show the validity of the 
above generalization. More recently, a direct verification of the Kadanoff's real space 
renormalization group picture of the Ising model in the case of a Sierpinski fractal has 
been achieved by Hsiao and Monceau |H1 ; they used a Monte Carlo renormalization group 
method. 

The simulation works mentioned above have been mainly performed with the help 
of the Monte Carlo Wolff algorithm 9]. This algorithm is able to reduce efficiently the 
critical slowing down. Instead of involving a single-spin flip at each update of the spin 
state configuration, the Wolff algorithm grows firstly a cluster (so-called Wolff cluster) 
and then updates the state configuration by assigning a new state to the spins of this 
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cluster. Much information on the critical behavior of a discrete spin system can be 
brought out from the geometrical properties of the Wolff clusters. It has been shown that 
the mean size of Wolff clusters for the Ising model scales as L'^/'^ at where 7 and v 
are the critical exponents associated to the susceptibility and the correlation length |[9|. 
It provides an alternative way to calculate 7/1/. Monceau and Hsiao checked that the 
values of 7/1/ calculated from the Wolff mean cluster sizes are consistent with the ones 
calculated from the behavior of the maxima of the susceptibility with respect to L in 
the case of the Sierpinski fractals with Hausdorff dimensions betweeen 1 and 3. [Tm ITT] 
Moreover, the most striking result they carried out is the scaling invariance of the Wolff 
cluster size probability densities V{s) under a suitable rescahng. This scaling law holds 
at the critical point and reads 

V{s) = 6-^^T(r"s), (2) 

where s denotes the size of Wolff cluster associated with a Monte Carlo step and h is 
some appropriate change of length scale which keeps invariant the fractal structure. They 
proved that the exponent x is equal to yh = /?/z/ + 7/z/ where (3 is the critical exponent 
associated to the magnetic moment. This scaling law can be understood as a decrease 
of Wolff cluster size by a factor h^^ under the change of the length scale 1^6. During 
the simulations of the Ising model, the total magnetization M and the total energy E 
are calculated at each update of the spin state configuration. As a matter of fact, the 
absolute difference | AM | between two successive values of M is equal to two times the size 
of the Wolff cluster grown during a Monte Carlo step. Hence, the scaling law described 
in Eq. (j2)) is satisfied by the probability density distribution of |AM|. It is worth noticing 
that the mean absolute difference |Aii^| between two successive values of the total energy 
represents two times the surface energy of the Wolff cluster. Thus, there is much to 
learn in studying the scaling properties of |Ai?| and the associated probability density 
distribution. The purpose of this paper is to study the behavior of the surface energy 
of Wolff clusters grown on fractal lattices at the critical point. We investigate seven 
different Sierpinski-type fractals of Hausdorff dimension lying between 1 and 3. The 
lattices are generated iteratively from some generating cells and denoted SP{£'^,nocc, k). 
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^ is the size of the generating cell, d is the integer space dimension in which the lattice is 
embedded, riocc is the number of occupied sites in the generating cell, and k is the number 
of iteration steps. Ising spins are placed at each center of the occupied sites. The size of 
SP{1'^, riocc, k) is L = and the number of spins is = n^^^,. The Hausdorff dimension is 
defined hj df = log A^/ log L = log nocc/ log ^. The true mathematical fractal is obtained 
only when k tends to infinity; in this case, we denote it by SP{i'^,nocc)- The structure 
of the fractal is not indicated in these symbols. The seven generating cells SP{£'^, Uocc, 1) 
we have chosen are 

(1) . ^P(22,4, 1): a 2 by 2 square, 

(2) . S'P(3^, 8, 1): a 3 by 3 square where the center sub-square is removed, 

(3) . S'P(5^, 24, 1): a 5 by 5 square where the center sub-square is removed, 

(4) . 5P(23,8,1): a cube of size 2, 

(5) . SP{3^, 26, 1): a cube of size 3 where the center sub-cube is removed, 

(6) . S'P(4^, 56, 1): a cube of size 4 where the center sub-cubes of size 2 are removed, 

(7) . 5'P(3^, 18, 1): one removes in addition the 8 sub-cubes at the corners of S'P(3'^, 26, 1). 

5'P(2^, 4) and SP{2^, 8) are exactly a square and a cubic lattice of infinite size. It has been 
shown that the Ising model exhibits a second order ferromagnetic phase transition on a 
fractal, provided that the lattice has a particular geometrical property: the ramification 
order must be infinite [12 . It is the case for the fractals studied here. Moreover, the 
critical temperatures Tc on these fractals are available [H 13 their values are 

recalled in the table 1. The mean absolute value of the Wolff cluster surface energy can 
be calculated from the relation 

where A^sim is the total number of Monte Carlo steps and En is the total energy of the 
n-th updated configuration. Firstly, we found that (|Ai?|) follows power laws at with 
respect to the lattice size L in the case of the seven different fractals we investigated. 



It enables to define a surface exponent: {\AE\) ~ L^^ . Fig.l sliows tlie beliavior of 
(lAi^l) as a function of L in logarithmic coordinates. The points line up along straight 
lines except for the small sizes, where the scaling corrections due to the finite-size effects 
are expected. It has been suggested that these corrections have a topological character 
and are linked to the slow convergence towards the thermodynamical limit in the case of 
the fractals with broken translational symmetry [Sj. As a matter of fact, (|A£'|) follows 
perfectly a power law in the case of the translationally invariant lattices SP{2^, 4, k) and 
SP{2^, 8, k) where L increasing as a geometrical series covers many orders of magnitude. 
We report the measured surface exponent Sw in the table 1, where least-square fits are 
performed from = 8 to 12 for S'P(2^,4, fc), from = 5 to 8 for SP{2^,8,k), from 
A; = 4 to 8 for SP{3'^,8,k), from A; = 2 to 4 for SP{4:\56,k), from A; = 3 to 5 for 
5P(52,24,A;), for SP{3^, 18, k), and for SP(3^26,A:), respectively. On the other hand, 
since Monte Carlo simulations can be performed only on lattices of finite size and we 
have omitted the small-size data when extrapolating the thermodynamical limit, a slow 
crossover behavior may be interpreted as an asymptotic one. A detailed study shows that 
the points in Fig.l exhibit a very slight concavity. The reported value Sw, hence, should 
be taken as an upper bound in a strict sense; however, the real thermodynamical-limit 
value is expected to be very close to the reported one. 

We are now able to go further and study the scahng properties of the probability den- 
sity distributions V{\AE\) of Wolff cluster surface energy. The curves showing V{\AE\) 
are similar to each other in logarithmic coordinates. These results suggest to write down 
a homogeneous transformation under the form: 

V{\AE\) = b-^^V{b-ys\AE\), (4) 

where Ds and ys are some introduced exponents. An intuitive guess is to set Ds equal 
to df — 1 since it describes the usual surface dimension of a system of bulk dimension df. 
According to Eq.(@]), {\AE\) should scale as : 

POO 

{\AE\) = / \AE\V{\AE\) d\AE\ 
Jo 

roo 

= / \AE\b'^'^f-^'^Vib-y'^AE\)d\AE\ 
Jo 

POO 

= h-^'^f-^^+^y^ / \AE'\V{\AE'\)d\AE'l 
Jo 
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where we have performed the change of variable \AE'\ = b~'^^\AE\. Since the correlation 
length is divergent at T^, we can set b equal to the lattice size L = (lAi^l) is, 
therefore, proportional to L~^'^f~^^^'^y'=^ . We are, hence, able to obtain a relation linking 
the exponents Sw and ys- 

Sw = -idf-l) + 2ys. (5) 

The values of ys calculated from Eq.(jni) for the seven fractals investigated are given in 
the table 1. 

The similarity property of the curves showing the probability density distributions 
V{\AE\) and the validity of the homogeneous transformation Eq.Q can be brought out 
in the following way: for a given structure at different values of the iteration step k, the 
curves showing V{\AE\) collapse onto the one corresponding to the lattice of the largest 
size L = under the mapping: 

Fig. 2 shows such collapses for the three largest values of k on the fractals S'P(3^,8), 
5'P(2^,4), and 5'P(3'^,18). These data-collapses work out in a reliable way with the 
values of ys given in Table 1. It confirms that Ds is equal to df — 1 (It has also been 
checked for the four other fractals). Furthermore, V{\AE\) does not exhibit a peak 
as V{s) does. The effect of segregation between large and small clusters mentioned in 
Ref. |im lTT] is smoothed in the behavior of V{\AE\). It means that an important part of 
the simulation is carried out in updating large clusters; however, the update of the large 
clusters does not necessarily imply a large change of the total energy. For a given cluster 
size, the probability distribution of the border surface is broad. 

A fundamental question is to know if the introduced exponent ys is a new one, that 
is, if it is independent of the two renormalization group eigen-exponents yt, yt and the 
Hausdorff dimension df. To our knowledge, it's the first time that the surface energy 
scaling property of Wolff cluster is investigated. No analytical or theoretical treatment 
is available. According to our Monte Carlo simulation results (table 1), ys seems to be 
independent of yt, yt and df. Whether ys links to some of the surface exponents of a bulk 
system [l5j remains an open question. In the table 1, regardless of the fractal structures, 
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one can find that the values of yt, Uh and ys increase as the Hausdorff dimension dj 
increases. Moreover, we do not expect that e-expansion results can be interpolated to 
non-integer dimensions and provide the values of yt and on the fractals, even the value 
of 2/5. Accurate Monte Carlo studies have shown that the universality of phase transitions 
on the hierarchical lattices without translational symmetry should depend on the lattice 
structure 01 Hence, ys associated to the fractals with broken translational symmetry 
depends on the lattice structure too. 

As the Wolff clusters are dynamical objects, we exploit here the connection between 
the surface exponent ys and the dynamical scaling exponent. We firstly study the mean 
square surface energy of Wolff cluster 

{mf) = j-^ E {E^-^i-Enf. (6) 

sim J- n=l 

We find that ((Ai?)^) lines up along straight lines with respect to L except for the small 
lattice sizes in logarithmic coordinates for the seven fractals investigated (see Fig. 3). 
The associated exponent m, defined by {{AE)"^) ~ L^", could be measured in the same 
way as described in the previous text for the exponent Sw', the values of u are re- 
ported in the table 1. We, then, express ((AE)^) as 2 [{E^) - (E)^] (1 - 9ee{1)) where 
(^ee{i^) = ((-Eo-E'n) (-2'^) ~ (-2')^) is the normalized autocorrelation function. Since 

the autocorrelation time is much longer than 1 (in some appropriate time unit), 1 — 6ee{^) 
represents approximatively the negative derivative of the autocorrelation function at the 
origin. It enables to define a statistic-fluctuation autocorrelation time and an asso- 
ciated dynamical exponent zfj-: —{d6EE{t)/dt\t=o) = i'^^f)^^ ~ L^^^f Notice that, 
according to Eq.(IT]), the term {E'^) — (E)"^ should asjTiiptotically scale as L'^/'^^^f = LF'/^ 
at Tc, where a is the critical exponent associated to the specific heat per spin c^,. One 
should keep in mind that the non-singular part of the free energy per spin gives, in ad- 
dition, an important contribution in determination of the critical behavior of 4j. We, 
therefore, have 2u = 2v^^ — zfj in the thermodynamical limit. With the help of the 
inequality (|Ai?|) < ((AE')^), we can link the surface exponent Sw with the dynamical 
one zfjr in the following way: Sw is upper bounded by m = u^^ — zfj-/2. Moreover, since 
{{AEY) follows a power law, a homogeneous transformation for the probability density 
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of {AEY can be stated; it takes the following form 

r{{AE)^) = b~^'^''f-^^V{b~y{AEf), (7) 

and has been verified on the seven fractals investigated here. In this case, {AEY decreases 
by a factor = under a suitable change of length unit 1 ^ 6. A similar process 

could be applied in the study of Wolff cluster size or the successive change of total 
magnetic moments. One can show that (|AM|) < ^ {[AMY) ^^"^5 therefore, 7/1/ is 
upper bounded by {jji' + df — Zgjr)/2 where zfj is the statistic-fluctuation dynamical 
exponent associated to the total magnetic moment. It implies z^jr < df — ^ / v = 2(3 / v. 

The method developed in this paper can be generalized to the study of some physical 
quantity whose mean value follows a power law with respect to L. For instance, since 
(|M|) ~ L^'' at Tc, one can verify, by constructing a homogeneous transformation for 
the probability density of \M\ similar to Eq.Q, that \M\ decreases by a factor b^y^+'^f)/'^ 
under a suitable rescaling 1 — > 6. All these results suggest that the scaling properties of 
some physical quantity at Tc, where the correlation length is divergent, comes originally 
from the underlying hierarchal structure of a fractal. 

In summary, we have studied the scaling properties of the Wolff cluster surface energy 
in the framework of the Ising model in the case of seven different fractals dimensions. A 
new scaling relation for the absolute value of the surface energy \AE\ of the Wolff's cluster 
has been established. We have shown that \AE\ scales as b^y^\AE\ under an increment 
of the length unit by the factor 6, which remains invariant the underlying structure of the 
fractal. Finally, the surface exponent Sw is proved to be upper bounded by — zfj/2. 
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Table 



Table 1: Measured values of Sw and u on the seven fractals, ys is calculated from 
Eq.©. The values of Tc, yt and yh are recalled from Ref. [S], Ref.jH], Ref.jTl], Ref.jlj, and 
Ref.gn]. 



fractal 




yt 


Vh 


Sw 


ys 


u 


SP(32,8) 


1.4795(5) 


0.449(6) 


1.8198(11) 


0.838(2) 


0.865(1) 


0.878(4) 


SP(52,24) 


2.0660(15) 


0.923(2) 


1.861(5) 


0.8188(2) 


0.8967(1) 


0.880(2) 


5P(22,4) 


2/ln(l + v^) 


1 


1.875 


0.8100(2) 


0.9050(1) 


0.8777(8) 


5P(33,18) 


2.35090(9) 


1.185(16) 


2.317(16) 


0.849(8) 


1.240(4) 


1.016(9) 


5P(43,56) 


3.99893(10) 


1.410(36) 


2.407(14) 


0.742(10) 


1.323(5) 


0.990(11) 


SP(3^26) 


4.21701(6) 


1.503(53) 


2.449(22) 


0.738(2) 


1.352(1) 


0.995(5) 


5P(23,8) 


4.511516(41) 


1.588(2) 


2.482(2) 


0.759(5) 


1.3795(25) 


1.018(5) 
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Figure captions 

Fig.l (|A£'|) versus lattice size L on the seven fractals investigated in logarithmic coor- 
dinates. 

Fig.2 Collapses of V{\l^E\) on the fractals 5P(32,8,A;), SP{2'^,A,k) and 5P(3^ 18, A;) 
under the mapping described in the text. 




lattice size L on the seven fractals investigated in logarithmic 



coordinates. 
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